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Abstract 

We study the three-body charmed baryonic decays B° — > T,^ + pir~ and i?° — > S[!p7r + in 
the four-body final state 5° — > A+p7r + 7r~, using a data sample of 357 fb -1 accumulated at 
the T(45) resonance with the Belle detector at the KEKB asymmetric-energy e + e~ collider. 
We measure the branching fractions B(B° -► £ c (2455) ++ p7r-)=(2.1 ± 0.2 ± 0.3 ± 0.5) xlO" 4 , 
B(B° -► S c (2455) V + )=(l-4 ± 0.2 ± 0.2 ± 0.4) x 10" 4 and B(B° -► E c (2520) ++ p7r _ )=(1.2 ± 0.1 ± 
0.2 ± 0.3)xl0~ 4 with signal significances of 13.1a, 9.4 a and 7.1cj, respectively. The errors are 
statistical, systematic, and due to the uncertainty in B(A^ — > pK~n + ), respectively. We also set 
an upper limit B(B° -» S c (2520)°p7r+) < 0.38 x 10~ 4 at the 90% confidence level. In addition, we 
obtain a non-resonant branching fraction of (6.4 ± 0.4 ± 0.9 ± 1.7) xlO -4 , and a total branching 
fraction of (11.2 ± 0.5 ± 1.4 ± 2.9) xlO" 4 for B° -» A+pn+Tr'. 

PACS numbers: 13.25.Hw, 14.20.Lq, 14.40.Nd 
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The large mass of the 6-quark enables B mesons to decay into two baryons with additional 
pions. Since the CKM matrix element \V c b\ [l[ is substantially larger than \V u b\, these 
baryonic decays preferentially proceed through b — > c transitions and produce final states rich 
in charmed baryons. CLEO pioneered the study of these processes and reported branching 
fractions and evidence of several exclusive charmed baryonic decays with a 9.1 fb _1 data 
sample @, 0, 0] • Recently, Belle has observed several new decay modes into two-, three-, and 
four-body final states with charmed baryons 0, 0], and three-body decays with charmless 
baryons p], 0, 0, EH HI- We find a hierarchy of the branching fractions that depends on 
the multiplicity in the final state: ~ 2 x 10 -5 for the two-body decays B° — > A+p and 
Br -> £°(2455/2520)p, ~ 1 x 10" 4 for the three-body decay Br -> A+pvr-, and ~ 7 x 10" 



for the four-body decay B° — > A^pn + n^ [6j, |l2j. There are several theoretical models that 
describe the decay mechanisms and predict the branchin g fr actions of baryonic B decays 
into two-body and three-body final states [H, 14, 15, Ila1l7l 18, 19, 20]. Detailed studies 
of such decays are very important to provide strict constraints on these theoretical models. 

In this paper, we report improved measurements of the intermediate three-body decays 
B° -> S c (2455/2520) ++ p7T- and B° -> S c (2455/2520) pvr + in the four-body final state 
B° — > A^p7i + n^. This study is based on a 357fb _1 data sample accumulated at the T(4S") 



with the Belle detector at the KEKB asymmetric-energy e + e~ collider [21 . 

The Belle detector is a large-solid- angle spectrometer based on a 1.5 Tesla superconduct- 
ing solenoid magnet. It consists of a silicon vertex detector (a three- layer silicon vertex 
detector (SVDI) for the first sample of (152.0 ±1.2) x 10 6 BB events and a four-layer silicon 
vertex detector (SVDII) for the latter (235.8 ±3.6) x 10 6 BB events), a 50-layer central drift 
chamber (CDC), an array of aerogel threshold Cherenkov counters (ACC), a barrel-like ar- 
rangement of time-of-flight scintillation counters (TOF), and an electromagnetic calorimeter 
comprised of Csl (Tl) crystals located inside the superconducting solenoid coil. An iron flux 
return located outside the coil is instrumented to detect K\ mesons and to identify muons. 
The detector is described in detail elsewhere (2~2| . 

We simulate the detector response and estimate the efficiency for signal reconstruction 
by Monte Carlo simulation (MC). We use the QQ program [23[ for signal event generation 
and a GEANT-based detector simulation program [24|. A sample of 5.45 x 10 4 signal 
events is generated for each of the four-body decay B° — ■> A^p7T + 7T~ , the intermediate three- 
body decays B° -> £ c (2455/2520) + Vr- and B° -»■ S c (2455/2520)°p7r + , and their charge 
conjugate modes. Each signal sample is processed by the detector simulation program that 
takes into account the differences between SVDI and SVDII as well as the long-term variation 
of the beam background conditions. 

The mode B° — > A^p7i^n + is tagged by an associated A+ particle, which decays into 
pK~n + . Charge-conjugate modes are implicitly included throughout this paper unless noted 
otherwise. To reconstruct A+ and B° signals, we require tracks to have distances-of-closest- 
approach to the interaction point of less than 5.0 cm in z (the direction opposite to the 
e + beam direction) and 1.0 cm in a plane perpendicular to the z-axis. We require the A+ 
mass to be within ±0.014 GeV/c 2 (~ 3.5cr) of our fitted mass of 2.287 GeV/c 2 . Hadrons 
such as protons, kaons and pions are identified by using likelihood ratios provided from the 
CDC dE/dx, TOF and ACC information (PID) @. We use likelihood ratios L S /(L S + 
Lf,), where s and b stand for the hadron species to be identified and for the background, 
respectively. We require these ratios to be greater than 0.6, 0.6 and 0.4 for proton, kaon and 
pion selection, respectively. The efficiency for proton identification is 95% with a kaon fake 
rate of 1.0%. The efficiencies for kaons and pions are about 90%; the corresponding pion 
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and kaon misidentification rates are about 10% 25j. Tracks that are positively identified as 
electrons or muons are rejected. We impose loose requirements on the vertex fit x 2 's for the 
tracks from A+ — > pK~ir + (x\+) an d B° — > A+p7r + 7r~ (Xb) to reject background from the 
decay products of K s and A particles. When there are multiple B candidates in an event, 
we choose the candidate with the smallest Xb- 

To suppress the continuum background (it, d, s, c pair production), we use two event shape 
variables, R 2 and cos thrust- The variable R 2 is the ratio of the second to zeroth order Fox- 
Wolfram moments [26[ , and cos thrust is defined as the cosine of the angle between the thrust 
axis of the reconstructed B decay products and the thrust axis of the other tracks in the 
center-of-mass system(CMS). We require R 2 < 0.35 and | cos thrust | < 0.8, which retain 80% 
of the signal and remove 60% of the continuum background. These requirements, together 
with the final B signal selection discussed below, reduce the continuum background by a 
factor of about 10 4 . 

The final selection requirements are based on the kinematic variables AE and M^ c . The 
variable AE = E B — E heSiin is the difference between the reconstructed B meson energy 



(Eb) and the beam energy (E^ eam ) evaluated in the CMS, while Mb c = J E^ eam — P| is the 
beam energy constrained B meson mass. Here, P B is the momentum of the B meson also 
evaluated in the CMS. 
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FIG. 1: (a) Scatter plot of AE vs. M bc for B° -> A+pyr+vr" signal candidates, (b) The AE 
distribution for M bc > 5.27GeV/c 2 . The shaded regions indicate the sideband whose total area is 
equal to the background in the B signal region, (c) The Mb c distribution for \AE\ < 0.03 GeV. 



Figured] (a) shows a scatter plot of AE vs. Mt, c . The vertical lines show the B signal 
region of 5.27GeV/c 2 < M^ c < 5.29GeV/c 2 , and the horizontal lines indicate the signal 
region, \AE\ < 0.03 GeV. Figure [T](b) shows the AE distribution for the M^ c signal region, 
where the curve shows the result of the fit with a double Gaussian for the signal and a 
linear background in the fit interval of —0.1 GeV< AE < 0.2 GeV. Figure [1] (c) is the M^ c 
distribution for |A.E| < 0.03 GeV. The curve shows the fit with a single Gaussian for the 
signal and an ARGUS function 27j] for the background. We use the AE distribution to 
determine the B signal yield, as we find a peaking background in the Mb c distribution from 
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a study of the AE sideband. 

To remove the feed-down from higher multiplicity modes with additional pions, we restrict 
the fit region to AE > — O.lGeV. The signal shape parameters are fixed to those fitted 
to the corresponding MC, where we find the Gaussian widths o\ (with a ratio of 02/01) of 
6.6±0.4 MeV/c 2 (2.2±0.2) for SVDI, and 7.2±0.2 MeV/c 2 (2.3±0.2) for SVDII. We obtain B 
signals of 535 ±30 and 865 ±38 events for SVDI and SVDII data, respectively. The efficiency- 
corrected signal yields normalized to the number of the BB events are consistent within 
errors. Thus, we combine the SVDI and SVDII data and obtain the total B° — > A+p7r + 7r~ 
yield of 1400 ± 49 events. 
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FIG. 2: The mass distributions of (a) A^7r + and (b) A^7r _ for B° — ► A^p7r + 7r _ . The points 
with error bars show the mass distribution for the events in the B signal region, and the shaded 
histogram indicates that for the sideband region. See the text for details. 



Figure [2] shows (a) the A+7r + and (b) the A+7r~ mass distributions for B° — ► A^p7i + n^. 
We observe clear peaks for the S c (2455) ++ '° and S c (2520) ++ . The points with error bars 
show the events in the B signal region defined by \AE\ < 0.03 GeV and 5.27GeV/c 2 < 
Mbc < 5.29GeV/c 2 . The shaded histograms are the events in the AE sideband region 
defined by 0.07 GeV < \AE\ < 0.10 GeV and 5.27GeV/c 2 < M bc < 5.29GeV/c 2 . As seen 
in Figure [1] (b), the number of the sideband events is equal to that of the background in 
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TABLE I: The Breit-Wigner parameters of £ c (2455/2520) resonances. In the simultaneous fits, 
they are fixed to Mpdg and rp^Q. Mg t and are the fitted values to the data with statistical 
errors only. 

Resonances M fit (MeV/c 2 ) fjf (MeV/c 2 ) M PDG (MeV/c 2 ) Tf DG {MeV/c 2 ) 



£ c (2455)++ 


2454.1 ± 0.2 


3.5 ±0.5 


2454.0 ±0.2 


3.44 ± 0.30 


£ c (2455)° 


2453.4 ± 0.6 


3.4 ±0.4 


2453.8 ± 0.2 


3.44 ± 0.40 


S c (2520)++ 


2517.9 ± 1.4 


19.9 ±3.5 


2518.4 ±0.6 


16.0 ±2.0 


£ c (2520)° 


2514.3 ±2.8 


19.1 ± 5.7 


2518.0 ±0.5 


16.0 ±2.0 



the B signal region of |A£?| < 0.03 GeV. We find small peaking backgrounds near the 
£ c (2455/2520)° masses in the sideband events. 

To obtain the £ c (2455) ++ /° and E c (2520) ++ /° signal yields, we consider possible contri- 
butions from peaking backgrounds seen in the sideband. We also study the background 
shape using MC samples for the non-resonant four-body decay B° — > A+p7r + 7r~, and the 
intermediate three-body decays B° -> S c (2455/2520) ++ pvr- and B° -> S c (2455/2520)°p7r + - 
We find a linear behavior for the mass distributions from non-resonant A+7T combinations. 
Therefore, we introduce independent linear background functions in the B signal and the 
sideband regions. 

We perform a simultaneous binned likelihood fit to the mass distributions with the fol- 
lowing functions 

(N ls + N lb ) x BWi(M) + (N 2s + N 2b ) x BW 2 (M) + (c s + a s x M) (1) 

N lb x BW 1 (M) + N 2b x BW 2 {M) + (c b + a b x M) (2) 

for the B signal and the sideband events, respectively. Here, M is the A+7T mass, and 
BW(M) represents a Breit-Wigner function. The subscripts 1 and 2 indicate S c (2455) and 
S c (2520), and the subscripts s and b stand for the signal and sideband, respectively, c and a 
are parameters of the linear functions. N\ s and N 2s are the net signal yields of S c (2455) and 
S c (2520), respectively, and N\ b and N 2b are the normalizations of the peaking backgrounds 
in the sidebands. 

Table fl] lists the parameters M PDG and Tp^Q for the Breit-Wigner functions BW\ and 
BW 2 used in the simultaneous fit. When we float these parameters, the fitter obtains values 
Met and T^, consistent with the PDG values M PDG and TpQ G 2^ (the latter is r PDG 



convoluted with the Belle detector resolution). Thus, we fix those parameters to the PDG 
values; the fitted parameters in the fit are then the signal yields A^,, and N 2s , the peaking 
background yields A^ 16 and N 2b , and the linear background shape parameters c s , a s , c b and 
a b . The uncertainties in the signal yields due to the assumed masses and widths are taken 
into account in systematic errors as discussed below. 

In the fits shown in Figure [21 we obtain x 2 /n.d.f = 183.4/192 and 196.6/192 for the 
fits to A+7r + and A+7r~ mass distributions, respectively. The solid curves show the fits to 
the mass distributions for the B signal region, and the dashed curves indicate the fits for 
the sideband region. The significance of the E c (2455) (S c (2520)) signal is evaluated as S = 



2 ln(L /L max ), where L max is the maximum likelihood of the fit and L is the likelihood 
for a fit with the yield of S c (2455) (S c (2520)) fixed to zero and the other parameters floated. 
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We study the change in the signal significances by varying the fixed masses and widths by 
their ±lcr errors and find that the resulting change is negligibly small. 

Table HU summarizes the fitted signal yields, efficiencies, significances and the branching 
fractions obtained for intermediate three-body decays B° — > E c (2455) ++< -°- ) £>7r~( + ) and B° — > 
£ c (2520)++( ( Vr~ (+) - The third error is due to the uncertainty in the branching fraction of 
£>(A+ — > pK~n + ) = (5.0 ± 1.3)%. As a check, we calculate separate branching fractions for 
charge-conjugate modes; the two branching fractions are in good agreement. 

We obtain a systematic error of 11.7% as a quadratic sum of 7.2% due to track recon- 
struction efficiency, 9.1% from the PID (both are coherent sums over the six tracks for the 
B decay products) and 1.9% due to the uncertainty on N(BB) and limited MC statistics. 
These errors are common to all decay modes. The signal efficiencies in Table [Til include the 
MC PID correction factor of 0.867 ± 0.079, to account for a systematic difference between 
data and MC. Separate PID correction factors for proton, kaon and pion tracks as functions 
of momentum and azimuthal angle are determined from a comparison of data and MC for 
large samples of D* + — > D°(K7i)'n- + and A — > decays. The overall PID correction factor 
is then calculated as a coherent sum over the six tracks for the selected B signal events. 
The error of ±0.079 is taken into account as the PID systematic error of 9.1% as mentioned 
above. We estimate an error of 3.5% for the total B signal yield from the maximum variation 
of the yield in fits to the AE distribution with the double Gaussian fixed to MC and with the 
shape paramters floated. This uncertainty in the B signal yield results in an error of 5.3% 
for the signal yield of the non- resonant four-body decay (see below). We estimate an error 
of 4.8% (9.1%) for S c (2455)++/° (S c (2520)++/°) from the variation in the fitted signal yield 
due to a ±lcr change (0.4 (2.0) MeV/c 2 ) in the width rpo G . We find a negligibly small effect 
on the mass. In addition, we take into account the uncertainty in the signal efficiency due 
to differences between the resonant substructure in data and signal MC. The S c (2455) ++ 
data is consistent with three-body MC, while the S c (2455)° data shows a broad p7r + mass 
structure that differs from MC phase space. We estimate an error of 4.6% for the E c (2455)° 
efficiency. 

We investigate the signal yield for non-resonant B° — > A+p7r + 7r~ decay, which consists 
of four-body decay, as well as contributions from decay modes with possible final state 
interactions or resonance states of the pir^, tt + tt~ and A+p systems. The signal efficiencies 
for p7r ± tend to be lower than that for the non-resonant four-body B decay near the mass 
threshold. We study two-body submass distributions and find some deviation from phase 
space near the threshold. However, due to limited statistics we cannot draw any strong 
conclusions about possible resonant structures. We conservatively estimate an uncertainty 
in the signal efficiency due to resonant structure to be 5%. Adding those errors in quadrature, 
we obtain total systematic errors of 12.6% for S c (2455) ++ , 13.5% for S c (2455)°, 14.8% for 
E c (2520)++/°, and 13.7% for non-resonant four-body B decay. 

We obtain the branching fraction of the non-resonant four-body decay by subtracting the 
signal yields for the observed three-body decays from the total B signal of 1400 ± 49 events 
and correcting for the efficiency of non-resonant four-body MC. The total branching fraction 
is obtained by adding the branching fractions of the intermediate three-body and non- 
resonant four-body decay modes. The branching fractions are consistent with the previous 
measurements (3, [1] and supersede our previous measurements (Bf . 

In summary, we study the three-body charmed baryonic decays B° — > E+ + p7r~ and B° — > 
YP c pir + in the four-body final state B° — > A+p7r + 7r~, and measure the branching fractions 
B(B° -> £ C (2455) + Vr-) = (2.1 ± 0.2 ± 0.3 ± 0.5)xl0~ 4 , B(B° -> £ c (2455)Vr + ) = (i-4 ± 



TABLE II: Branching fractions for B° -> £ c (2455/2520) ++ p7r~ and B° -» S c (2455/2520)°pvr + . 
The errors in the branching fractions are statistical, systematic and due to the uncertainty in 
B(Af — > pET _ 7r + ) = 5.0 ± 1.3%, respectively. See text for details of the systematic errors. 





Modes 


Yield 


Det.eff.(%) Sys.err(%) Si 


gn.(o-) 


^(xlO" 4 ) 


B°- 


-> S c (2455) ++ p7T- 


182 ± 15 


4.57 


12.6 


13.1 


2.1 ±0.2 ±0.3 ±0.5 


B° 


-» S c (2455) p7T+ 


122 ± 14 


4.41 


13.5 


9.4 


1.4 ±0.2 ±0.2 ±0.4 


B°- 


-> S c (2520) ++ p7T- 


155 ± 18 


6.91 


14.8 


7.1 


1.2 ±0.1 ±0.2 ±0.3 


B° 


-» S c (2520)°p7T + 


22 ± 16 


6.75 


14.8 


1.3 


< 0.38 (90%C.L.) 




Non-resonant 


919 ± 58 


7.50 


13.7 




6.4 ± 0.4 ± 0.9 ± 1.7 




Total 


1400 ± 49 








11. 2 ±0.5 ±1.4 ±2. 9 



0.2±0.2±0.4)xl0~ 4 and B(B° -> S c (2520) ++ p7r-) = (1.2 ±0.1 ±0.2 ±0.3) x 1(T 4 with signal 
significances of 13. la, 9.4a and 7.1a, respectively. The errors are statistical, systematic, 
and due to the uncertainty in £>(A+ — > pK~ir + ), respectively. We also set an upper limit 
B{B° -> S C (2520)V + ) < 0.38 x 10~ 4 at the 90% confidence level. In addition, we obtain 
a non-resonant branching fraction of (6.4 ± 0.4 ± 0.9 ± 1.7) xlO -4 , and a total branching 
fraction of (11.2 ± 0.5 ± 1.4 ± 2.9)xl0~ 4 _for B° -> A+pn+n-. 

The observed branching fraction B(B° — > S c (2455) ++ p7r _ ) is comparable to a previ- 
ous measurement of B{B~ — > A+p7v~) by Belle [29(. The S c (2455) ++ mode has a larger 
branching fraction than the E c (2455)° and S c (2520) ++ modes, and the S c (2520)° mode is 
significantly suppressed. The branching fraction B(B° — > S c (2455/2520) ++ pir~) is larger 
than B(B° — > S c (2455/2520)°p 7r+ ) ) probably due to an additional contribution from an 
external W emission diagram [20J. The total branching fraction for the four-body decay 
B(B° — ► A+p7T + 7r~) is five times larger than B(B~ — ► A+p7r~), as it consists of both inter- 
mediate three-body decays B° — > E c (2455/2520)p7r (~40%) and a non-resonant four-body 
decay (~60%). 
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